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, Abstract 

The process of the vacuum polarization energy losses of high energy cosmic rays 
propagating in the extragalactic space is considered. The process is due to the polar- 
ization of Cosmic Background Radiation by a moving charged particle. With the goal 
of the description of the process, the photon mass, refractive indices and permittivity 
■ function for low and high energy photons are found. Calculations show the rather 

noticeable level of the energy losses for propagating protons with the energies more 
\ than 10 — 10 GeV. The influence of the polarization energy losses on propagation 

7— I ■ of cosmic rays is discussed. 

\o : 
o . 

1 Introduction 

Q_i! The problems of the origin and propagation of high energy cosmic rays is widely discussed 
in recent years (see, for example JI[ 0, [| and literature therein). The photoproduction 
of mesons and e ± -pairs in the Cosmic Background Radiation was considered as the main 
reason of the energy losses of charged particles propagating in the extragalactic space 0, 
||. However, the simulations of the spectral distributions of cosmic rays on this basis 
do not provide a close agreement with observed data 0, and some peculiarities in the 
spectrum like the "knee" || and GZK-cutofT |7], || do not find recognized explanation. In 
this connection the determination of another reason of the energy losses is important and 
interesting. 

In this paper we consider the possible mechanism by which cosmic rays lose their 
energy in the extragalactic space. It is the polarization energy losses of a charged particle 
moving in the electromagnetic vacuum. In the presence of an external electromagnetic 
field the polarization of vacuum was considered first in the pioneer papers || . The vacuum 
polarization leads to different effects ]10|] such as the nonlinearity of the Maxwell equations, 



appearance of nonzero photon mass, birefringence of light, etc. The description of the 
various methods in the QED of vacuum one can find in literature ( see |TT| and references 
therein). In parallel with other methods the traditional method, based on the introduction 
of the permittivity tensor, may be used for the effective description of different phenomena 
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in the electromagnetic vacuum. In this case the vacuum described by the Maxwell equations 
which are similar to equations of the classical electrodynamics for continuous media JT2], |T5| . 
As an example, one can point out the well known low energy permittivity and permeability 
tensors for constant electric and magnetic fields J1U[|. In paper the connection between 



the polarization and permittivity tensors is established. The using of the permittivity 
tensor is convenient because this approach allows to consider the QED vacuum and matter 
with the unit point of view. 



2 The mass of photon, propagating in an isotropic 
photon gas 

The photon mass m 7 in an isotropic medium is defined by the refractive index n of 
monochromatic photons propagating in this medium 

my = 2£?(1 - n(E 7 )) (1) 

where E^ is the photon energy. From this equation one can see the equivalence of finding of 
the photon mass and the refractive index. Note that in the general case the refractive index 
is a complex value. However, in a transparent medium its imaginary part is a negligibly 
small quantity. 

Our calculations of photon mass in an isotropic photon gas are based on the results of 



papers |L5|, |16| , where propagation of 7-quanta in the field of monochromatic and dichro- 
matic laser waves was considered. In these papers the permittivity tensor was derived for 
such media. Knowing this tensor one can find the refractive indices of photon and other 
characteristics of the propagation process. The imaginary components of the tensor are 
connected with the energy losses due to e ± -pair production in the laser wave. The real 
components of the permittivity tensor are derived with the help of dispersion relations. 
In the case of an isotropic photon gas the energy losses one can describe by the following 
relation (see, for example |T7|]): 



m; a x 
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where L 77 is the reciprocal of the mean free path for collisions, n(e) is the differential 
photon gas number density for the photon energy equal to e, s is the center of momentum 
frame energy squared, cr 77 (s) is the total cross section of e ± -pair production, E^ is the 
energy of photon propagating in a medium, s min = Am 2 c A , s max = 4e£ 7 , e min = s min /4:E 7 , 
m is the electron mass and c is the speed of light in the vacuum. All relations in this paper 
will be valid for the laboratory coordinate system, which determined as the system where 
the mean momentum of the background photons is equal to zero. From Eq.(2) one can 
find the imaginary part of refractive index n 

MW = (3) 
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The explicit form of this relation is: 

nr 2 X 



Im(n(E y )) = — j J 7 im (z m )n(e)ede, (4) 



Fim(Zm) = Zl f [(1 + Z- Z 2 /2) 111 \ + V L^- ~ VT^(1 + Z )}^ (5) 

Jz m 1 - vl - z z A 

where r e and A e are the classical radius of electron and its Compton wavelength, z = 
i.m 2 c A /s and z m = 4m 2 c 4 / ' s max = m 2 c A / '(eE 7 ). The corresponding real part of the photon 
refractive index is: 



At a r°° 

Re(n(E 1 )) - 1 = e 9 e / J r re (z m )n(e)ede ) (6) 
mc 2 Jo 

where function T re {z) is determined by the following relation: 

T re {z m ) = zl i r ^{z)- 2 (7) 



-3 - (1 + z - ^)\L 2 _ -(l-(z+ z i))\L\+ 

+ ^£^L_ - + £ (1 + z - f ), < z < 1, 

-3 + (1 + z - 4) axccot 2 (v^T) - (1 - (z + ^ ))\L\+ 1 J 



k +(1 + ^v^Iarccot sfz~=\ - {z ~ 1) / W * L + , z > 1. 
The functions L_, L+ are equal to: 



L_=ln±±*EI, L + = ln^±±. (9) 

The obtained here relations for refractive indices are valid for low and high energy photons 
propagating in a photon gas. Now one can calculate the refractive indices of the photon 
propagating in the extragalactic space. As the first approximation of this medium one can 
take the model of the space which is filled by the Cosmic Background Radiation. In this 
case the number density is defined by well known equation: 

1 e 2 de 

n{e)de = 3 — , 10 

■K 2 c 6 h exp(e/T) - 1 

where T is the temperature. For low energy photons one can find the following relation 
for real part of the refractive index: 

Re(n(0)) - 1 = / n(e)ede = I ® \ * ■ (11) 

y y " 135mc 2 Jo v ; 2025 m 4 c 8 v ; 

The real and imaginary parts of refractive indices are connected by the dispersion relation: 
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At present temperature T = 2.726° K the real part of refractive index for low energy 
photons is 5.1 10~ 43 . These results are in agrement with the similar calculations [15 , |TP| , |2"0| 
of the low energy refractive indices in an isotropic photon gas. 

Fig.l illustrates the results of calculation of the photon refractive indices in the simplest 
model of the extragalactic medium (in the Cosmic Background Radiation). One can see 
that the real part is nearly constant in the wide energy range from till 2 10 5 GeV with 
the flat maximum (« 6 1CT 43 ) at m 3 10 5 GeV. At energies more than 2.5 10 6 GeV 
the real part of index is negative. The imaginary part has maximum value ~ 4 10~ 43 at 
£ 7 = 8 10 5 GeV. 

Knowing the refractive indices one can find the energy dependence of the permittivity 
function of a photon gas. However, this function is not uniquely defined and it depends 
on the relation between magnetic induction vector B and intensity of magnetic field H 



see |L2], [13]]). In |15[ the Maxwell equations are written with H = B and in this case the 



permittivity function e is 

e{E 1 ) = n 2 « 1 + 2Re(n - 1) + 2ilm(n). (13) 



3 Energy losses of particle moving in a medium 

Propagating in the Cosmic Background Radiation charged cosmic rays produce polarization 
of this medium. As result the cosmic rays lose the initial energy. This process is similar 
in many respects to the ionization losses of energy of charged particles moving in matter. 
It is clear that the vacuum polarization is many times more weak process, however it may 
be noticeable on a large distances of particle propagation in such a specific medium as the 
extragalactic space. 

Our consideration will be based on the Landau theory of ionization energy losses of 
relativistic particles in a matter |T2[. However, a small adaptation of the theory will be 



done with the goal of the elimination of divergences in the final result. 
The Maxwell equations in a medium have the following form: 

1 f)Yi 

divB = 0, rotE = — , (14) 

c ot 

I deE 4tt. 

div£E = 4vrp, rotB = — 1 j. (15) 

c ot c 

where E is the intensity of electric field and B is the magnetic induction vector, p and j 
are the charge and current densities, e is the permittivity operator (see|T^]), t is the time. 
The charge and current distributions one can take in the following form: 

p = ep(r — vf), j = ep(r — v£)v, / p(r)dr = 1, (16) 

Jv 

where e is the particle charge. The field potentials we write in the usual form: 

B = rotA, E = --— -grady? (17) 
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In accordance with [12] we use the further condition for potentials: 

divA+i^ = 0. (18) 
c at 

The substitution Eq(17) in Eq(15) allows to get the following relations for potentials: 

. , e d 2 A Air , , . 

AA "^^ = "T ep(r " vt)v ' (19) 

i(A V - = -4nep(r - vt). (20) 

These equations for Fourier components have the form: 

F f) 2 Ai 4-7T 

k 2 Ak + = — eF k vexpHkvt), (21) 

£ d 2 

e{k 2 ^ + 4^Jr) = 47reF k exp (-ikvt). (22) 
cr ot z 

where F k — J p(r) exp(— ikr)c?r. One can see that the functions A k and (p depend on the 
time as exp (— ikv). As it shown in [jn|, the result of action of e-operator on exp (—icut)- 
function is its multiplication on e(o;)-function. Now from Eqs. (21-22) one can get: 

Attc F^v 

Ak = T P - (kv)%(kv)/c 2 GXP ( " Zkvt) ' (23) 
Aire i^k 

^ = I(kv) P - (kv)%(kv)/c 2 6XP { ~ lkvt) - (24) 
The Fourier component of the intensity of electric field has the following form: 

zkv 

E k = A k - zk<p k . (25) 

c 

Now we find the damping force, which act on the extended charge. For this purpose we 
use the following well known relation: 

|(^) = -JE-divS, S = ^[EB], (26) 

where D is the electrical induction vector. After integration over some volume V one can 
get: 

^ = -ev / p(r)E(r)rfr - / S n da, H = f ° E Q + — dV, (27) 

dt JV Jo JV 07T 

where a is the bounding surface. It is clear, that at r — > oo the surface integral tends to 
zero. Then the damping force has the following form: 

T =-— [ p(r)E(r)rfr (28) 
v Jv 
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Obviously that Fouirier component of electric intensity has the form: = *fr(k) exp — ikvt, 
where <&(k) is the known function (see Eqs. (23-25)). Then we can write the following re- 
lation: 



ev / p(r- vt)Erfr = (29) 
Jv 

ev f r ev r 

J {y^p(r-vt)expzk(r-vt)rfr}$(k)rfk = ^-p J F k $(k) dk, 

Obviously, the da mping fo rce is directed against the velocity of particle. Let us name 
k x v = u> and q = <Jky + k 2 , where x-coordinate axis is along the direction of motion, to is 
the photon frequency. Now one can obtain the following equation for the absolute value of 
the damping force: 

T= — \ / k \ ^ 2( c l } \l 30 
7i J-^Jo e[q 2 + u} 2 {^ - js)] 

One can see that these relation is differ from similar relation in |IIJ by F 2 (/c)-multiplier. 
This multiplier is the charge distribution of the moving particle in k-space (the axial 
symmetry of the charge distribution around x-axis is assumed). For next we can write 
that e sa 1 + Ae' + iAe", \As'\ << 1, Ae" << 1. Then we extract the real and imaginary 
parts of the integrand function in Eq.(30). Thus, we get the sum of two integrals: 7i + Hi- 
Besides, it needs to take into account that Ae' and Ae" are even and odd functions of the 
frequency. One can show that I2 = 0, and therefore the damping force is: 

2e 2 /-oo ,oo F 2 (q 2 + -^)Ae"(uj)(q 2 + ^)ujqdujdq 



T ™ 2 Jo Jo (g 2 + 7 ^-^T) 2 + 7 %^ (31) 

where 7 is the Lorentz factor of the particle. Besides, in this equation F is the charge 
form factor of the particles (i.e. the charge distribution in k-space in the particle rest 
frame). The Eq.(31) for the damping force is final. The two integrals in this equation are 
converged. The integral over q is converged due to the form factor of the particle, and the 
integral over uj is converged due to relations — >• 1 at uj — >• 00. 



4 Energy losses of cosmic rays 

For cosmic rays Eq.(31) one can simplify. It is well known that the maximum observed 
energy of cosmic rays is less than 10 12 GeV. One can see that energy losses depend on the 
real and imaginary parts of the permittivity function. However, one can neglect of this 
dependence in the denominator. It is really that 7 -2 >> |Ae'|, Ae"(see Fig.l and |2"I]). 
Now one can write the damping force for cosmic rays in the following form: 

rtWo Jo if+&r 

where q = hq and = Two. Here we use relation: e" = Ae". For the next calculations we 
take the empirical electric proton form factor 



w) = (i + ||r 2 - (33) 
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where Q is the three-dimensional transfer momentum, the empirical constant M v = 
(0.84GeV7c) 2 . Then we find the integral over q and get the following relation: 

T = \ ^'(E^S^E,) + S 2 {E 1 ))E 1 dE 1 (34) 

7m v 2 Jo 

y3 y2 

Si(-Ey) = - ln(xy) - y - — - y, (35) 

S 2 (E 7 ) = x(\ - l)(41n(x2/) + ^ + y 2 + 3y + -), (36) 

i? 2 1 

'j 2 v 2 M v 1 + x 

It is helpful to obtain the inexact but simple estimation relation for the damping force. It 
is possible to make for the large 7-factors (> 10 s ) of the cosmic rays. This relation is: 

T ~ -t-i ln (7^) [ El ' b e'XE^E.dE, (38) 
Tin v 2 < Ey > Jo 

where < E 1 > is some value of the 7-quantum energy range, where the integrand has 
noticeable quantities, E 1)b is the boundary value of the energy range and it is satisfied the 
condition E^^ 3> jMyV. For example, one can define the E 1 by the following equation: 

rE lb poo 

£ 7 = e"(E y )E^dE y / J e"(E 7 )dE T (39) 

According to calculations < E 7 >~ 8 10 7 GeV and the quantity of the integral in the 
numerator of Eq.(39) is equal to ps 10~ 28 GeV 2 at E 7ife ~ 10 13 GeV. The choice of E ljb in 
the energy range from 10 10 till 10 14 GeV is weakly changed the result. Then we get the 
following estimation: 

aZ 2 9a , 'jMyC , , 

irnc < E 7 > 

where Z is the particle charge in the elementary one units. From here one can see that 
energy losses of high energy cosmic rays (protons) is about some hundreds of GeV per light 
year. This value is the same order as the e ± -pair production in the Cosmic Background 
Radiation. Fig. 2 illustrates the calculation of vacuum polarization energy losses of protons 
in accordance with Eq.(34). 



5 Discussion 

The results of calculations of the vacuum polarization energy losses of the high energy 
cosmic rays show that losses are reasonably large, and, because of this, it is necessary 
to take into account this process for true description of the particle propagation in the 
extragalactic space. Fig. 3 illustrates the relative polarization energy losses of protons and 
iron nuclei in the extragalactic space. The results of calculation of the energy losses 
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due to the photoproduction processes (p + 7 — *p + e + + e~,p + r y — * p + tt°, ...) are 
also shown on the figure. One can see the following peculiarities of the considered process: 

i) the spectral distribution of polarization energy losses are rather broad, and the flat 
maximum of the relative losses is at E pmax ps 4 10 7 GeV and E Fe>max rs 1.4 \Q w GeV for 
protons and iron nuclei, correspondingly; 

ii) the relative value of the polarization energy losses has the same order as the losses 
due to the photoproduction processes (at proton energies in the range 10 7 — 10 10 GeV); 

iii) the strong drop of polarization losses take place for particle energies <C E Ptmax , Ep e ,max; 

iv) the spectral behavior of proton and iron polarization losses is the same, but there 
is valuable energy shift between the curves on Figs. 2-3. As result at energies < 10 9 GeV 
the polarization losses of protons is more, than ones for iron nuclei. The reverse situation 
is observed at energies > 10 9 GeV; 

v) our calculations show that the points of inflection are exist for the both curves on 
Fig. 2 (when d?T/d 2 E = 0). It takes place when the energies of particles are equal to 
1.5 10 7 and 1.4 10 10 GeV for protons and iron nuclei, correspondingly. 

The behavior of energy losses for iron nuclei one can understand if taken into account 
that the losses are proportional to the square of charge but the form factor cuts of the value 
of damping force more strong for iron, than for proton (the constant My = 0.13GeV/c ). 

Now we make attempt to explain some experimental data in the observation of the 
cosmic ray spectrum. The absence of the clear photoabsorbtion threshold in the spectrum 
H |5| one can explain by the common continuous character of the summary energy losses. 
Really, from fig. 3 we see that the polarization energy losses sewed together with the pho- 
toabsorbtion ones. It means that clear photoproduction threshold can not be detected on 
this background. 

The " knee" -effect is another misunderstand area in the spectrum of the cosmic rays. 
One can suggest that the "knee" is a point where the character of cosmic ray propagation 
is changed. At energies above the "knee" -level the particles lose their relative energy and 
after passing the "knee" -point the energy losses are decreased and accumulation particles 
in this area takes place. 

Some remarks concerning the composition of cosmic rays. The composition is dom- 
inated by protons at the lowest energies, and then the fraction of light nuclei increases 
with energy ^3], |24|. However, at energies > 5 10 9 GeV the protons is again dominated. 
We can see on Fig. 2 that at low energies the polarization losses of protons in many times 
exceed the same losses for iron nuclei. On the other hand, at energies > 5 10 9 polarization 
losses for iron nuclei are large. In particular, the composition of cosmic rays is determined 
by lifespan in which a cosmic particle keeps the energy. Note, the "knee" is observed only 
for proton fraction, and it is absent for iron nuclei with energies < 10 9 GeV. From our 
point of view this fact is obvious. We believe that the iron "knee" is at energies more than 
5 10 9 GeV near the point of inflection for the iron energy losses curve. 

In the case if the considered here mechanism of vacuum polarization energy losses is 
true the following statement take place: the initial number (or the speed of production) of 
cosmic rays ( with energies above the "knee"-point) is more, than it is commonly supported. 

These our conclusions have the qualitative character. However, no doubt they may be 
tested by simulation of the propagation process of cosmic rays. 

In conclusion, we have touched on the nature of the considered phenomenon. The 
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using of methods of the electrodynamics of the continuous media allow to describe math- 
ematically the mechanism of the vacuum polarization energy losses of relativistic particles 
without detailed description of the primary processes, which are responsible for this phe- 
nomenon. From this standpoint the moving charged particles polarize the medium. It 
requires some energy and the particles give back its to the medium. In the case of the elec- 
tromagnetic vacuum this energy go into creation of virtual e ± -pairs. According to Eq.(39) 
the spectrum of the virtual pairs is the rather broad and its upper bound is near the energy 
of particles. Although the quantities of the permittivity function e" are rather small, the 
noticeable value of the energy losses is generated due to broad and high energy spectrum 
of virtual pairs in the laboratory coordinate system. It is obvious that these pairs are low 
energetic in the rest frame of the cosmic charged particle. Then the real photons from the 
Cosmic Background Radiation and virtual ones can interact effectively in between at the 
condition eE lv > rr^c 4 , where E 1V is the energy of the virtual state. 

It should be noted, that we do not consider the influence of the red shift on the energy 
losses. It needs to make in the case when propagation of the cosmic rays is investigated on 
long distances. Besides, we think that the contribution in the vacuum polarization energy 
losses of cosmic rays from another background fields is possible. 



6 Conclusion 

On the basic of determination such characteristics as the photon mass, refractive indices 
and permittivity function in the Cosmic Background Radiation the vacuum polarization 
losses of high energy cosmic rays are considered. The calculations show the high level of 
these losses for protons with energies more than « 10 7 GeV. The proposed mechanism 
of losses leads to a revision of the existing propagation models of cosmic rays. With our 
point of view the propagation of the high energy cosmic rays in the extragalactic space is 
the dynamic process to a greater extent than it is expected. Experimental and theoretical 
ivestigations of these processes will help to understand the nature and origin of the cosmic 
rays in the Universe. 

The author would like to thank H. Zaraket for critical questions, remarks and useful 
references. 
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Re(n)- 1 




Figure 1: The real minus unit (1) and imaginary (2) parts of refractive index (in 
linear (a) and logarithmic (b) scales) in Cosmic Background Radiation as functions 
of the photon energy. The absolute value \Reh — 1| is shown in logarithmic scale. 



11 



10 5 10 7 10 9 10 1 E.GeV 



Figure 2: The vacuum polarization energy losses of protons and iron nuclei propa- 
gating in the extragalactic space as function of the energy. 



12 



-7 

10 - 




Figure 3: The relative vacuum polarization energy losses of protons and iron 
nuclei propagating in the extragalactic space as function of the energy. On the 
right the similar value for photoproduction energy losses (from [2]) is presented for 
comparison. 
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